Electronic Theory for Bilayer-Effects in High-T c Superconductors 
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The normal and the superconducting state of two coupled CuOi layers in the High-T c super- 
conductors are investigated by using the bilayer Hubbard model, the FLEX approximation on the 
real frequency axis and the Eliashberg theory. We find that the planes are antiferromagnetically 
correlated which leads to a strongly enhanced shadow band formation. Furthermore, the inter-layer 
hopping is renormalized which causes a blocking of the quasi particle inter-plane transfer for low 
doping concentrations. Finally, the superconducting order parameter is found to have a d x 2_ y 2 
symmetry with significant additional inter-layer contributions. 
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The importance of the multiple CuOi layers within 
the High-T c superconductors like Bi2Sr2CaCu20s+5 
(BSCCO) or YBa 2 Cu 3 6+5 (YBCO) is intensively stud- 
ied. Neutron Scattering experiments show clear evi- 
dence that antiferromagnetically correlated spin fluctu- 
ations arc quite strong between nearest-neighbor layers 
in YBCO f[l[§. Recent angular resolved photoemission 
(ARPES) experiments found indications for two sepa- 
rated bands in YBCO that might be related to 
the existence of two CuOi bands caused by a coherent 
inter-plane quasi particle transfer. However, the small 
experimental observed bilayer splitting in YBCO and the 
difficulty to resolve two Cu0 2 bands in BSCCO §0] 
support the idea that the strong short ranged antifer- 
romagnetic order in the cuprates alters the electronic 
excitations and reduces the inter-layer hopping. Here 
the observation of shadows of the Fermi surface (FS) in 
BSCCO demonstrates impressively that precursors 
of the magnetic phase transition are already present in 
the normal state. Furthermore, it has been suggested 
that a reduced quasi particle hopping and a Josephson- 
like coherent Cooper pair tunneling might be responsible 
for the interesting layer-dependence of T c jl0| . 

Theoretically, two coupled CuOi layers have been 
studied microscopically by incorporating solely a mag- 
netic coupling [fLl) and within the bilayer Hubbard model 



Despite these efforts it remains an open question to de- 
termine the influence of the bilayer coupling on the quasi 
particle dispersions, on the inter-layer charge dynamics, 
on c-axis transport properties and most of all on the su- 
perconducting phase of the cuprates. Here, it might be 
of particular importance to investigate the relation be- 
tween electronic theories and the results of the model 
Hamiltonian of Chakravarty et al. jl0| . 

In this letter we study the bilayer Hubbard model 
within the FLEX approximation [Q on the real fre- 
quency axis. We find that the planes are antiferromag- 
netically correlated yielding strong deformations of the 
quasi particle dispersions and a blocking of the effective 
inter-layer hopping for small doping concentrations and 
low excitation energies. For YBCO- and BSCCO-like 
systems we observe that shadow states occur only when 
the inter-plane coupling is considered. For the supercon- 
ducting state we find that the d x i_ y i symmetry is most 
stable and that a bilayer is characterized by inter- and 
intra-layer Cooper pair formation. 

Our theory is based on the Hubbard Hamiltonian for 
a bilayer CwCVsystem 
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including a finite inter-plane hopping 
recent phenomenological theories jl4|, 15 1 , it has been ar- 
gued that the presence of bilayer correlations might lead 
to a superconducting s-wave order parameter that has 
opposite signs in the two C11O2 bands (s -state). More- 
over, for a single layer compound with a model dispersion 
we recently presented results indicating that the shadow 
states Q are caused by an antiferromagnctic coupling of 
momenta k at the FS and k + Q (Q = (vr, 7r)) jUi and 
studied their importance for the fine structure in tunnel- 
ing and ARPES experiments below T c [^7|. However, it 
is a priori not clear how the realistic dispersions of the 
cuprates and the significant bilayer correlations affect the 
electronic excitation spectrum in BSCCO and YBCO. 



E4JI3II. Here and in where the hopping integrals t- v determine the bare dis 



persion £»/(k) in 2D k-space, i and j (I and I') are 
the site (layer) indices, the Kronecker symbol, U 
the local Coulomb repulsion and fi the chemical poten- 
tial. The interaction-free contribution of the Hamiltonian 
can be diagonalized yielding an antibonding (— ) and a 
bonding band (+) with bare dispersion e±(k). Assum- 
ing that this symmetry holds also in the full interact- 
ing case, one can define corresponding electronic Greens 
function G±(h,iuj m ) with fcrmionic Matsubara frequen- 
cies Lo m = (2m + l)nT and temperature T. They are con- 
nected to the self energy £±(k, iw m ) via the Dyson equa- 
tion G±(k,iuj m ) = l/(iu> m - (e±(k) -/x) - E±(k, iw m ))- 
For the FLEX self energy one obtains 
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with ^2 k , = m , , band index X = ± and number 
of momentum points iV. Note that due to the symme- 
try of the bilayer it follows for the intra-band interac- 
tions that V++(k, iuj m ) = K-_(k, iw m ) and for the inter- 
band contributions that V-\ (k, iuj m ) — V_+(k, iuj m ). 

V\ } \i (k, ioj m ) can be obtained by calculating the electron- 
hole bubble XA,A'(k,iw m ) = - J2k' G * ( k ' + k > iuJ m + 
iuj m ') G\i (k', iu) m i) and by performing the straightfor- 
ward summation of the FLEX diagrams |Q in the layer 
representation yielding an inter-plane, Vj_(k, iu) m )i and 
an in-plane, Vii (k, iu) m ), contribution: 
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FIG. 1. Dependence of the effective interactions on the 
doping and on the inter-layer hopping t±. (a) Ratio -Vj_/Vii 
and V++/V+- for the LSCO-like model, x = 0.09 (solid 
line), x — 0.12 (dashed line) and x — 0.16 (dashed-dotted 
line). Fermi surfaces and their shadows of the LSCO- (b) and 
YBCO-like model (c) dispersions. 

This set of coupled equations is solved self-consistently 
on the real frequency axis by using the numerical method 
presented in Ref. |l9| . To characterize the FS of the dif- 
ferent High-T c superconductors we use the bare disper- 
sion 

£±(k) = — [2t(cos(fc x ) + cos(ky)) + At' cos(k x ) cos(k y ) 
+2<"(cos(2fc K ) + cos(2fc y )) ± tjj 

with hopping integrals t — 0.25 eV, t' = —0.38 t and 
t" = -0.06 t as model for YBCO and compare 

these results with a model dispersion for if = t" = 0. 
Since the later FS is similar to the La2_ 2; Sr 2 ;Cu04 sys- 
tem for t± = 0, we call this dispersion for simplicity 
LSCO-like although LSCO has in distinction to YBCO 
or BSCCO only one plane per unite cell. We investigate 
the complete dependence of the electronic properties on 
t±, but take t± =QAt= 100 meV as inter-plane hopping 



in YBCO [g3|]. For the comparison with previous results 
we take U = 4t but notice that we find no significant 
changes in our data up to values of U = 6t. 

In Fig. 1(a) we investigate the interplay between 
the inter- and intra-plane antiferromagnetic correlations 
in bilayer cuprates upon t± and the doping x, where 
n = 1 — x is the occupation number per site. By 
using the LSCO-like dispersion we find that Vji (k, us) 
decreases only slightly when t± is increased and that 
the magnetic in-plane correlation length £ f[& is al- 
most not sensitive to the inter-layer hopping p2j . Since 
Vj_(k, u>) and Vji(k, w) have the same k-dependence, but 
with opposite signs, the two layers are antiferromagnet- 
ically correlated as observed in experiments As a 
measure of this phenomenon we plot in Fig. 1(a) the 
ratio of the maxima of the corresponding interactions 
— Vj_/Vj| := — Vj_(Q, u))/Vu(Q, uj). Here we observe that 
the bilayer coupling increases rapidly with increasing t± 
until it saturates. In addition, it is remarkable that not 
only Vj_(k, to) and Vji (k, w) increase with decreasing dop- 
ing, but also the ratio — Vj_/Vji. Thus, the stabilization 
of the inter-plane magnetism via a coupling of antiferro- 
magnetically ordered in-plane regions (size ~ £) across 
the layers should be most effective for low doping con- 
centrations. By transforming these results to the band 
representation via Eq. [l], one sees in Fig. 1(a), where we 
plot V++/V+- := V++(Q,w)/V+_(Q,w) ) that the con- 
tribution of the intra-band interactions vanishes rapidly 
for low doping with increasing t± thereby explaining the 
absence of these excitations in experiments [Q . The sup- 
pression of V++(Q,ui) and the corresponding enhance- 
ment of K| (Q,w) can also be understood by consider- 
ing the FS topology of the LSCO- and YBCO-like models 
in Fig. 1(b) and (c), where the bonding and the anti- 
bonding FS are almost shifted by Q and the fact that 
we find very flat quasi particle bands at (tt, 0). Finally, 
the later point leads also to an antibonding FS that is 
closed around (0,0), which is in contrast to LDA calcu- 
lations p3|] and might also be interesting for the current 
debate about the FS in YBCO §. 

The dominance of the inter-layer and inter-band corre- 
lations has important consequences for the quasi particle 
dispersions and the bilayer splitting. By investigating the 
bonding and the antibonding bands of the LSCO model 
with t± = 0.4 t [^U, we observe similar to Ref. |l3} | 
that the bilayer splitting is strongly renormalized from 
the uncorrelated Ae(k = (tt,0)) = 2 tj_ = 200 meV to 
Ae(k = (tt,0)) = 125 meV in the case of x = 0.12. To 
study the origin of this effect and more interestingly to 
investigate the low energy inter-planc electron dynamics 
Re S_L(k,w) = 0.5 Re (£+(k,cj) - E_(k,o;)) is plot- 
ted in Fig. 2 for k = (tt, 0) upon doping. This quan- 
tity renormalizes tj_ and yields an effective frequency 
and momentum dependent inter-plane hopping ampli- 
tude ij_(k, u) = t± — Re £_i_(k, w), which follows directly 
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from the matrix Dyson equation in the layer representa- 
tion. The sign of Re £j_ (k, to) is determined for small to 
by the sign of Vj_(k, to), such that we find ij_(k, to) <C t± 
for to < t± and low x (x < 0.12). This interesting phe- 
nomenon can be found along the entire FS although it 
is most pronounced at (n, 0) and exists for a wide range 
of t± as can be observed in the inset of Fig. 2. Con- 
sequently, the effective inter-plane hopping at the Fermi 
energy is blocked in the cuprates, which leads to the dis- 
appearance of a coherent quasi particle tunneling process 
although due to the frequency dependence of Re S j_ (k, to) 
both bands are still splitted. Physically, this is related to 
the fact that each inter-layer hopping process is accompa- 
nied by spin-flips due to the opposite antiferromagnetic 
environment in the other plane. Thus, due to the cor- 
responding large Im Sj_(k, to) an incoherent coupling of 
the layers occur. Note that a hypothetical ferromagnetic 
coupling across the bilayers would not suppress tj_(k, to), 
since there is no competition between the kinetic energy 
gain and the magnetic correlations. Finally, the remark- 
able doping dependence of our results might be of par- 
ticular importance for the c-axis transport properties of 
the cuprates, since underdosed and overdosed compounds 
are characterized by qualitatively different temperature 
dependences p5|. 
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FIG. 2. Re E x (k,cj) for k = (tt, 0) and for different doping 
values. Inset: Re £±(k, 0) at the FS of the antibonding band 
versus t±, where uj — refers to the Fermi energy. 



In Fig. 3 we demonstrate the influence of the inter- 
plane coupling on the formation of shadow states. In 
Rcf. we showed for a single layer system that the an- 
tiferromagnetic coupling of momenta k and k + Q leads 
to a transfer of spectral weight from the FS to its shadow. 
Now, since we find for a bilayer system that V++ (k, to) 

is small compared to V-\ (k, to) for intermediate t±, we 

expect that the spectral weight is not only shifted by the 
momentum Q, but simultaneously also from the bond- 



ing to the antibonding band and vise versa. This phe- 
nomenon is demonstrated in Fig. 3 for the LSCO-like 
model dispersion, where we present our results for the 
spectral density for the bonding (g + (kps, to)) and the 
antibonding (g-(kps + Q, to)) band and FS momentum 
k^g. There is a transfer from the bonding to the anti- 
bonding band leading to an occupied shadow state and 
in addition a transfer in the reversed direction above the 
Fermi energy. Most interestingly, the intensity of the 
shadow states is strongly dependent on the magnitude of 
the inter-plane hopping as can be observed in the right 
inset of Fig. 3. 
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FIG. 3. Spectral density Q±(k.FS,u) for the LSCO-like 
model. Right inset: Dependence of the shadow state intensity 
on t±. Left inset: Position of the bands at the (7r,0) point 
with respect to the Fermi energy for (7 = and U ^ 0. 

Here the shadow peak is most pronounced for t± = 
0.4 t before it decreases again for even larger values of t±. 
To verify that this result is related to an enhanced antifer- 
romagnetic coupling, we investigated Im £_(kp g + Q, 0) 
as a measure of the coupling of the FS to its shadow |lq | 
and found a corresponding maximum at tj_ = 0.4 t. 
Therefore, the strongly enhanced formation of shadow 
states in bilayers is caused by electronic correlations and 
is related to the already discussed suppression of the bi- 
layer splitting and a corresponding large spectral density 
at the Fermi level. This is demonstrated in the left inset 
of Fig. 3, where we plot the rjj_-dependence of the posi- 
tion (w( ff) o)) of the flat quasi particle band at (7r,0) with 
respect to the Fermi level for U ^ and U = 0. It can 
be seen that the antibonding band moves away from the 
Fermi energy (~ t±), while the bonding band is tightly 
pinned at to — up to t± « 0.45 t. Note that the shadow 
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state intensity is also enhanced for larger x and most im- 
portantly that in a bilayer shadow states can be found up 
to doping concentrations of x rj 0.17, while they vanish 
for x f» 0.13 in a single layer compound. 

By performing calculations for the YBCO dispersion 
we find no shadow states for t± = 0. Therefore, due to 
the similarity of the YBCO and BSCCO FS |§7| the ex- 
perimental observation by Aebi et al. || can not be sat- 
isfactorily understood by considering only a single CuOi 
plane. However, by taking a finite t± into account we find 
for YBCO in agreement with our results for the LSCO- 
like model that the antiferromagnctic coupling increases 
and that shadow states start to appear for x < 0.12 with 
a maximum intensity at t± w 0.4 t. Furthermore we 
find that the most favorable region to observe shadow 
states in ARPES is in the neighborhood of the (tt/2, tt/2) 
point, where main and shadow band are well separated. 
Near (jr, 0) the absolute intensity of the shadow states 
is largest, but they are difficult to detect because of the 
superposition of shadow peaks and the dominant main 
band contributions ]27| ]. 

(a) d-wave, single layer (b) s^wave, bilayer 




(c) d-wave, intra layer (d) d-wave, inter layer 
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FIG. 4. Superconducting order parameter for uj — 0, 
T = 48 K and U = At (in meV). (a) single layer (t ± = 0). 
(b) meta-stable bilayer s -state for the bonding band and 
t± = 0.4 t. (c-d) bilayer d-wave state in the layer representa- 
tion for t± = 0.4 t. 

The superconducting state of the bilayer Hubbard 
model was treated by using a strong coupling Eliash- 
berg theory within the fully self-consistent FLEX approx- 
imation ]2q| . By assuming that there is no inter-band 
pairing, one obtains one order parameter for each band, 
namely <fi± (k, ui) , which are connected to the layer repre- 
sentation via 4>± (k, ui) = (k, ui) ± 0j_ (k, ui) . By solving 
these equations for the YBCO dispersion with t± = 0, we 
find for all x a superconducting state with d x 2_ y 2 pairing 
symmetry and values of T c (T c = 64 K for x = 0.12) that 
are only slightly smaller than the T c of the previously in- 
vestigated LSCO-like dispersion (T c = 78 K) 0. Con- 
cerning the pairing symmetry in a bilayer system, we ex- 
pect in view of our results for the effective interactions for 



larger doping an intra-plane d x 2_ y 2 state. However, since 
the inter-band interactions are dominant for low doping 
concentrations, the argumentation of Ref. O] might ap- 
ply leading to an s^-symmetry, which corresponds to an 
inter-plane pairing state. To clarify this point, we present 
in Fig. 4 our results for T well below T c and the LSCO- 
like model with tj_ = 0.4 t and x = 0.09 in comparison 
with t± — |^2j . Here in contrast to the single layer case, 
where only a d-wave pairing state occurs (Fig. 4(a)), we 
find an s ± -state as a solution of the Eliashberg equations 
as shown in Fig. 4(b). Note, this order parameter is 
highly anisotropic and largest at the FS to achieve the 
highest energy gain. Nevertheless, the s ± -state refers 
only to a meta-stable solution, since the bilayer d x i_ y i 
pairing symmetry, as shown in Fig. 4(c) and (d), still 
yields cf> d 3> 4> s and consequently a much larger conden- 
sation energy. Interestingly, the <i-wave state is charac- 
terized by an increasing contribution of inter-layer pair- 
ing with decreasing doping, where for x = 0.09 Cooper 
pairs are formed by electrons from the same layer as from 
different layers with almost equal probability. However, 
we did not find a large enhancement of the correspond- 
ing T c , which might be caused by magnetic frustrations, 
since the inter-layer d-wave pairing occurs also for lattice 
sites that are ferromagnetically coupled above T c . 

In conclusion, we presented new results for the dop- 
ing dependence of the bilayer Hubbard model by using 
the FLEX approximation on the real frequency axis. For 
realistic inter-plane hopping integrals we find that the 
antifcrromagnetic coupling across the layers is strong in 
YBCO, which leads to an enhancement of the shadow 
state intensity and to the appearance of shadow states 
in YBCO- and BSCCO-like systems. Furthermore, the 
short ranged antiferromagnetic order reduces the bilayer 
splitting and blocks the inter-layer hopping for low dop- 
ing and small excitation energies. In the superconduct- 
ing state we found an order parameter with d x 2_ y 2 sym- 
metry that has equally important inter- and intra-layer 
pairing contributions. However, it remains an interest- 
ing problem to analyze the observed increase of T c with 
the number of layers within one unit cell. This might be 
achieved by including intrinsic Josephson tunneling into 
our description of the inter-layer d-wave state |]l0| . 
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